Onsager-Machlup theory and work fluctuation theorem for a harmonically driven 

Brownian particle. 
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We extend Tooru-Cohen analysis for nonequilirium steady state(NSS) of a Brownian particle to 
nonequilibrium oscillatory state (NOS) of Brownian particle by considering time dependent external 
drive protocol. We consider an unbounded charged Brownian particle in the presence of oscillating 
electric field and prove work fluctuation theorem, which is valid for any initial distribution and at all 
times. For harmonically bounded and constantly dragged Brownian particle considered by Tooru and 
Cohen, work fluctuation theorem is valid for any initial condition(also NSS), but only in large time 
limit. We use Onsager-Machlup Lagrangian with a constraint to obtain frequency dependent work 
distribution function, and describe entropy production rate and properties of dissipation functions 
for the present system using Onsager-Machlup functional. 
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I. INTRODUCTION 

The first theory of fluctuations in equilibrium thermo- 
dynamics was given by A. Einstein, who used the Boltz- 
mann relationship between the probability of a state and 
it's entropy. It states that the probability of a fluctu- 
ation from equilibrium state is proportional to e As ^ k ' B , 
where AS is the variation of entropy calculated along a 
reversible transformation creating the fluctuation and kg 
is the Boltzmann constant. 

In 1953, L. Onsager and S. Machlup generalized the 
equilibrium Boltzmann-Einstcin formula to near equi- 
librium dynamical setting (near equilibrium irreversible 
processes) by using a variational principle to calculate 
the most probable trajectory followed by the system in 
the spontaneous emergence of a fluctuation[l|,[f 
The laws of correlations in average course of a fluctu- 
ation (a sequence of statistically correlated nonequilib- 
rium states) are precisely given by the laws of irreversible 
processes. Under the assumption of time reversibility of 
microscopic dynamics, they showed that most probable 
creation and relaxation trajectories of a fluctuation are 
time reversal of one another. So, the path integral con- 
cept or the Onsager-Machlup functional gives a natural 
generalization to the time domain of Einstein-Boltzmann 
formula. 

Recently, fluctuation theorems has drawn considerable 
attention in nonequilibrium statistical mechanics. Fluc- 
tuation theorems are the statements about the asym- 
metry of the distribution functions of work, heat etc. 
around zero, i.e., they involve negative tails of work 
or heat. These relations points towards rare events in 
macroscopic systems, which are more easily observable 
in small systems with finite degrees of freedom. These 
theorems give more complete picture of the second law 
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of thermodynamics L 

known that near equilibrium fluctuation-dissipation re- 
lation and Onsager reciprocal relations can be obtained 
from fluctuation theorems [Tol [l9| . 

In a recent work[38j Tooru Taniguchi and E. G. D. 
Cohen extended Onsager-Machlup theory of fluctuations 
around equilibrium to a theory of fluctuations around 
nonequilibrium steady state for a specific system, con- 
sisting of a constantly dragged Brownian particle (B- 
particle) through fluid. They discuss thermodynamics 
of this system using generalized Onsager-Machlup La- 
grangian. They also discuss fluctuation theorems based 
on their generalized Onsager-Machlup theory and gen- 
eralized forms of the detailed balance conditions for 
nonequilibrium steady states. A very general treatment 
of macroscopic fluctuation theory for a class of continu- 
ous time Markov chains is given by Bertini etal[12, HH], 
and a general treatment considerin g la rge deviation the- 
ories is given by B. Derrida etal}34l. |35||. 

The present paper has a limited objective, in this 
work we extend Tooru-Cohen analysis for nonequilirium 
steady state(NSS) of B-particle to nonequilibrium os- 
cillatory state (NOS) of B-particle by considering time 
dependent external drive protocol. We consider an un- 
bounded charged Brownian particle in the presence of 
oscillating electric field (figure 1) and prove work fluctu- 
ation theorem, which is valid for any initial distribution 
and at all times. We use path integral approach similar 
to Tooru-Cohen approach [38j|. We use Onsager-Machlup 
Lagrangian with a constraint to obtain frequency depen- 
dent work distribution function. 

The paper is organized as follows. In section II we in- 
troduce our model and finite time transition probability 
using path integral approach. We then calculate it's most 
probable path using variational principle and solve func- 
tional integral for finite time transition probability. In 
section III we calculate work distribution function and 
prove work fluctuation theorem. In section IV we use 
Onsager-Machlup functional to describe coarsed-grained 
B-particle thermodynamics. In section V we end with 
conclusion. 
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CHARGED BROWNIAN PARTICLE IN 
OSCILLATING FIELD 




FIG. 1: Charged Brownian particle in an oscillating electric 
field 



II. A CHARGED BROWNIAN PARTICLE IN 
AN OSCILLATING ELECTRIC FIELD 

Consider a charged Brownian particle in an oscillating 
electric field as shown in figure 1. We assume that the 
external harmonic drive time scale is much larger than 
internal dynamical time scale ( inverse molecular collision 
frequency) of the B-particlc-fhrid system, i.e., 



Transition probability for a finite time interval [see Ap- 
pendix B for a complete calculation] is given by 



F 



Xi 

t 



J'O 

to 



Dx s exp 



dsL{x s ) 



(7) 



The Lagrangian L for this driven Brownian motion is; 

1 



L(x s ) 



AD 



(x s + 7] cos tos) , 



and the Path Integral is introduced as; 
ft 



Dx s exp 



dsL(x s ) 



it, 



1 



lim 



(8) 



N 



< J dx tN _ 1 J dx tN _ 2 ...J dx tl exp J 



dsL(x s ) 



(9) 



Here the initial time is to, the final time is t^ — t, and 
the initial position is x Q , the final position is x t , with 
t n = t + nAt, n= 1,2,3, ...,N, At = (t - t )/N. 



Solution Path: Calculation of the above 
Functional Integral 



> 



uj molecular collision frequency 



(1) 



Thus friction coefficient is independent from drive 
frequency [40l| . The B-particle's driven stochastic dynam- 
ics is given by the following Langevin equation, 



d 2 x t 
1 dt 2 



-a—— - qE cosujt + & (2) 
dt 



Here, £t is a Gaussian- White Noise(GWN) with proper- 
ties (&) = and = 2a/f3S(t - t'). In an over- 
damped case, it is 



dx t 
dt 



qE 



COSOjt 



if,. 

a 



(3) 



The Fokker-Planck Equation for the above Langevin 
equation is [see Appendix A for derivation] 



df d 

- = - VCOS Uf 



d 2 f - 



(4) 



Here / is the noise averaged distribution function, D = 
-^p is the diffusion coefficient, with j3 — l/fc^T, and r\ - 
gEa _ Ec_ _ rpjjg Fokker-Planck operator C is defined as 



fid r,d 2 
C = —ncosujt + D—r. 

ox ox 1 



(5) 



The transition probability for B-particle for a small step 
from (x' , t) to (x, t + At) is given by 



F 



X 


x 1 ' 


t + At 


t 



= e CAt S(x - x') 



(6) 



In this stochastic dynamical setting where the Brown- 
ian particle is fluctuating about the most probable path 
governed by external harmonic drive, the variational 
properties of the action governs the most probable tra- 
jectory of the particle between two given fixed points of 
the path. In the present setting the variational princi- 
ple takes care of two aspects of the B-particle motion, 
(1) the systematic part due to harmonic drive and (2) 
the stochastic part due to molecular impacts (which fur- 
ther involves fluctuation and dissipation). Is is to be 
noted here that the dissipative aspect of the motion 
which originates from two different physical effects (1) 
due to stochastic noise(fluctuation-dissipation), and (2) 
due to external systematic drive, are taken into account 
by phenomenological parameter a in the Langevin equa- 
tion, and in the present treatment we are not taking into 
account the frictional force which is nonlinear in velocity 
possible for an arbitrary drive strength. The variation 
principle which maximizes the transition probability i.e., 
for the solution path({i s }, from xq to Xt) is 



5 [ dsL{x s ) = 0, 

Jto 



which implies Eulcr-Lagrange equation 

d_ [9L\ dL 
ds \dx s J dx 



= 0. 



(10) 



(11) 



The solution path obtained from above the equation with 
end points at xq — xq and at x t = %t is given by, 



V ■ 

x s = smcos 



x t + ^ smujt — Xq 

t 



Xq. 



(12) 
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To calculate the functional integral in equation (9), we 



will first calculate the integral 



consider a variation z s 

■ .s 



J dsL(x s ) 



For this, 



x a about the solution path 



dL 



1 



,d 2 L 



dsL(x s + z s ) = I ds[L + z s ^- + ^(z s ) 2 —^ 
t Jto dx s 2 dxj 

(13) 

Inserting various derivatives from the known Lagrangian 
and solution path, the above integral is 

If 1 7] 

dsL(x s ) / ds[2(x t H — sin ujt — x ) + z s ]z s . 

The functional integral (equation 9) then becomes 



Xt 
t 



Xq 
to 



/ Dz s exp 

'zo 



1 

"Id 



ds(2g{t)/t + z s )z s 



to 



(15) 



Here g(t) = Xt + ^j smut— xq. Next we calculate the func- 
tional integral in the above equation(Eq.l5) for which 
zo(to = 0) = z t = 0. Writing the exponential part as 



exp 



~ f ds({i s +g{t)/tf-{g{t)/tf) 

J to 



exp 



1 



N 



ADAt 



- Zi + &tg(t)/t)f 



e 4Dt 



(16) 



Ef =1 (5 l+ i-2.) 2 -4fe(^+i)At(4^) 2 ] e ^r. (i 7 ) 



On taking the limit At — > or N — > oo with i — to = 
TVAi, the first part of the exponential after path integra- 
tion gives standard result l/vwi\Di(for standard solution 



see|41j). and the second gives e 



we finally obtain 
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Xt 


X() 


t 


to _ 



1 



^jAitDt 



e iD 



jyj (xt— sinwt) 2 



(18) 



Which reduces to free B-particle motion for rj = 
0(no external drive). 



III. FUNCTIONAL INTEGRAL CALCULATION 
FOR WORK DISTRIBUTION FUNCTION 

The rate of work done by the external oscillating elec- 
tric field is given by 



W{x t ) = (—qE cosu>t)xf 



(19) 



For simplicity we consider the dimensionless work 
(3Wt({x s }) and it's distribution defined as 

Pw(W,t) = {{8{W - (3W t ({x s })))) t . (20) 



Here, ((. . .)) t means a functional average over all possible 
paths {x s } between starting and end points as well as 
integrals over all starting and end points of the path: 

P w {W,t) = fdxt I * Dx s f dx f{x a ,t a )e ! ° dsL(± ° ) 



(21) 



x[S(W - f3W t ({x s }))}. 



The term e$° dsL i Xs ) is the probability functional for the 
path which emanate from any initial distribution denoted 
as f(xo,to). Inserting integral representation for delta 
function in the above equation (Eq. 21) we get 

P w {W,t) = — / d\ dx t Dx s / dxo 

xf(xo,to)eti d8L ^e iXW - iX P w *« x «» (22) 

Inserting the expression Wt({x s }) = f Q dsW(x 8 ) for 
work done in a finite time interval in the above equa- 
tion, we obtain 

Pw(W,t) = 
i r + °° 

— / d\e* xw 
2tt 



dx t I dx f(xo,to)F(x t ,x ;iX). 

(23) 



Where F(xt, Xq; A) is the constrained transition proba- 
bility 



F(x t ,xo;X)= / Dx s exp\ / ds[L(x s ) - Xf3W(x s )} 

Jxq VJo 

(24) 

With, L{x s ) — XpW(x s ) = L total- Next we calculate 
the solution path for the total Lagrangian from Euler- 
Lagrange equation, and then calculate the constrained 
transition probability F(xt, Xo\ A). Solving the Euler- 
Lagrange equation (Eq. 11) for the above total La- 
grangian we obtain the solution path as 

77(1 - 2A) . (x t + v(1 ~ 2X) 8mujt-x \ 
x s = — - smws+ ^ s+x . 

V 1 J 

(25) 

Proceeding on exactly similar lines for the calculation of 
transition probability (from equation 13 to 18) we obtain 
constrained (taking work rate into account) transition 
probability as 

F(x t ,x ; A) = 

^exp{-^(. 9l (t,A)) 2 -A, 2 (t) + A^ W 

, 7 ? „• 



g\(t, A) = Xt — xo + (1 — 2A)— sinajt 
, . n 2 , sin2cji. 



(26) 



Which reduces to the previous result (equation 18) for 
A = 0. The calculation for Pyy(W, t) is straightforward, 
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inserting equation 26 into equation 23, and using stan- 
dard Gaussian integrals and normalization condition for 
initial distribution, we have 



P w (W,t) = 



^£(1 + ^)* 



(27) 

From which we immediately obtain work fluctuation the- 
orem: 



Pw(W,t) 
Pw(~W,t) 



(28) 



The important point is that we have proved the work 
fluctuation theorem for any initial distribution and for all 
times as an identity (for a free driven B-particle, where 
the electrostatic force acts directly on the particle [3^]). 
For the system (a harmonically bounded and constantly 
dragged Brownian particle) considered by Tooru and 
Cohen[38], work fluctuation theorem is valid for any ini- 
tial condition(also NSS), but only in large time limit. 



IV. ONSAGER-MACHLUP FUNCTIONAL AND 
COARSE-GRAINED BROWNIAN PARTICLE 
THERMODYNAMICS 

From equations (7) and (8), the Onsager-Machlup 
(OM) Lagrangian function is defined as 



and both are positive functions. For noise averaged ve- 
locity (from noise averaging of Langevin equation) both 
the dissipation functions are equal to each other, 

1 F 2 

*((is) noise) = HF e , to) = — S- cos 2 LUS > 0. (33) 

21 a 

From Eq. (32) and Eq. (33) we obtain 

S((x) noise ) = 2$((i s ) noise ) = 2<&(F e , u) > 0. (34) 

In other words the sum of dissipation functions is equal 
to the rate of total entropy production. 

Next we will show that the energy conservation law 
is true for any fluctuation (without coarse-graning). The 
rate of work done on the B-particle by the external field 
is defined as W(x s ) = — F e (x s ) cos uis, and the total 
work done for a time interval t is given by PFt({a; s }) = 
J Q dsW(x s ). The heat generated Q t ({x s }) for the solu- 
tion path {x s } is defined as 

Qt({x s }) =T dsS(x s ) = — I dsXf,Fp cos ws 



(35) 



dsW(x s ) = W t ({x s }). 



Thus for a steady oscillatory state, the work done by the 
external field is dissipated as heat into the fluid surround- 
ing the Brownian particle. 



L(x s ) 



■ cos ojs 



(29) 



where F e = qEo. The above OM Lagrangian can be 
expressed as 



L(x s ) 



where S(x s 



1 

2fc^ 



*(x a ) + 9(F B ,u)-S(x.) 



(30) 



■k{xs)F e cos ujs is the Entropy produc- 



tion rate, <&(± s ) = Yr(x s ) 2 and Hf(F e ,Lj) 



IT 



COS LOS 

are Dissipation functions. The entropy production rate 
S(x s ) is a positive function for a coarsed grained descrip- 
tion 

Coarse - graning = (. . .) noise . (31) 
Using Langevin equation (Eq. (3)) we obtain 



F 2 

S((x) noise ) = cos 2 cut > 0, 



(32) 



which justifies the property of entropy production rate. 
It is clear from the definition of Dissipative Functions 

$(i s ) = ^(x s ) 2 and ^(F e ,u>) = ^y-^cos 2 cos that they 
are invariant under the time reversal change x s — > — x s , 



V. CONCLUSION 

For an unbounded charged Brownian particle in the 
presence of oscillating electric field we have proved work 
fluctuation theorem, which is valid for any initial dis- 
tribution and at all times. We have used path inte- 
gral approach similar to Tooru-Cohen approach and ex- 
tend their analysis for noncquilirium steady state(NSS) 
of a Brownian particle to nonequilibrium oscillatory state 
(NOS) of Brownian particle by considering time depen- 
dent external drive protocol. For harmonically bounded 
and constantly dragged Brownian particle considered by 
Tooru and Cohen, work fluctuation theorem is valid for 
any initial condition, but only in large time limit. We 
used Onsager-Machlup Lagrangian with a constraint to 
obtain frequency dependent work distribution function, 
and used Onsager-Muchlup functional to describe ther- 
modynamic properties of the present system. A complete 
treatment of this problem valid at any distance from equi- 
librium(when frictional force is nonlinear in velocity) and 
at all frequencies [10] is an open issue. 
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Appendix A: Fokker-Planck Equation 

To drive the Fokker-Planck Equation 42] for the 
Langevin equation(Eq.(3)), we define probability distri- 
bution of variable x at time t as f(x, t). Probability con- 
servation gives J dxf(x,t) — 1, or 



dt dx 1 Bt ' 



(36) 



| = -^(-,cos^/ + ie t /) 07) 

Define Lf = -§^{— rycos u)t)f, with rj — First we 

will consider the noise free part of the above equation 
i.e., ^ = — Lf, it's solution is 

/(*,*) = e"J> £ /M). (38) 
With noise term, 

%— t '-hb>- m 

we have a solution of the form 

f(x, t) = e~ /o dt ' L f(x, 0) - /* dse- IS dt ' L e K dt " L 

Jo 

x~S„f(x,8) (40) 
ox a 

Substituting Eq. (40) back into Eq. (39), and taking 
average over noise, we obtain 



d(f(x,t)) 
dt 

J_d_ '■' 

a 2 dx 



---L(f(x,t)) 

dse -iti dt'i-ss dt"i] d s)) 

ox 



(41) 



It is important to note that f(x,t) depends on noise £ s 
only for times s that are earlier that t. The RHS of 
the above equation contains two explicit noise factors, 
£t and £ s , and also those earlier time noise factors that 
are implicit in f(x, s). The first pair gives S(t — s), and 
second pair with the implicit noise factor gives S(t — s'), 



which does not contribute to the above integral as t > 
s > s' . Only first pair will contribute. Using properties 
of Gaussian white noise and delta function, we obtain 



df 
dt 



dx 



d 2 f 

rjcosujtf + d -q^2- 



(42) 



Here / is the noise averaged distribution function and 
diffusion coefficient D = l/a/3, with (3 — 1/fcsT. 

Appendix B: Transition Probability 

The transition probability for B-particle for "a small 
time step" At from (x' , t) to (x,t + At) is given bv[43| 



F 




X 


x' 


t- 


f At 


t 


X 


X 1 ' 






t + At 


t 


= cxp 



5{x — X ) 



(43) 



F 

xS(x — x') 



d . ^ A d 2 

^—At-q cos Lut + DAt—^r 

dx dx 2 



i\(x—x ) 



1 f +oc d d 2 

= — / dAcxpf— At-ncosujt + DAt— -z\e 
27r J_ 00 dx dx z 

-| /--(-oo , 
= — rfAe [--DA 2 At+ 4 AtA(£^+^cosu;t)]^ /^s 

27T — 



Standard Gaussian integral for A gives 

1 

F 



X 


x'~ 


t + At 


t 



y/AnDAt 



(45) 

Using Chapman-Kolmogorov equation the transition 
probability for a finite time interval is given by 





Xq 






t 

ft 


to 






Dx s exp 




/'■' * 










Dx s 


exp 


' J-'u 








) = 


1 

iD 


(x 





1 



t 

dsL(x 



ds(x s + 7] cosws)^ 



(46) 



Here, L(x s ) = —-^{x s + t^cosws) 2 is the Lagrangian for 
this driven Brownian motion. 
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